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In this work we show that the gravity lagrangian f(R) at relatively low curvatures in both metric 
and Palatini formalisms is a bounded function that can only depart from the linearity within the 
limits defined by well known functions. We obtain those functions by analyzing a set of inequalities 
that any f(R) theory must satisfy in order to be compatible with laboratory and solar system 
observational constraints. This result implies that the recently suggested f(R) gravity theories with 
nonlinear terms that dominate at low curvatures are incompatible with observations and, therefore, 
cannot represent a valid mechanism to justify the cosmic speed-up. 
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Observations carried out in the last few years, indicate 
that the universe is undergoing a period of accelerated ex- 
pansion 0. In the context of this cosmic speed-up, mod- 
ified theories of gravity in which the gravity lagrangian 
is a nonlinear function of the scalar curvature 



We find that the lagrangian must be linear in R and 
that the possible nonlinear corrections are bounded by 
R 2 . This implies that the cosmic speed-up cannot be 
due to unexpected gravitational effects at low cosmic 
curvatures. 
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have become object of recent interest. These theories are 
commonly referred to as f(R) gravities. Motivated by 
the fact that the addition of positive powers of the scalar 
curvature to the Hilbert-Einstein lagrangian may give 
rise to early-time inflationary periods, modifications of 
the lagrangian that become dominant at low curvatures 
have been suggested to justify the observed late-time 
cosmic acceleration (see also @). The aim of these 
theories is to describe the cosmic acceleration as an effect 
of the gravitational dynamics itself rather than as due 
to the existence of sources of dark energy. In addition 
to the standard formulation of f(R) gravities, where the 
metric is the only gravitational field, it was pointed out 
in Q that once nonlinear terms are introduced in the 
lagrangian, an alternative and inequivalent formulation 
of these theories is possible. If one considers metric 
and connection as independent fields, i.e., that the 
connection in the gravity lagrangian is not the usual 
Levi-Civita connection but is to be determined by 
the equations of motion, then the resulting theory is 
different from that defined in terms of the metric only. 
This new approach, known as Palatini formalism, and 
the metric one have been shown to lead to late-time 
cosmic acceleration for many different f(R) lagrangians. 
However, the attempts made so far to unravel the form 
of the function f(R) from the cosmological data are far 
from being conclusive 

In this work we study the constraints on f(R) gravities 
imposed by laboratory and solar system experiments. 



In order to confront the predictions of a given gravity 
theory with experiment in the solar system, it is nec- 
essary to obtain its weak field, slow motion (or post- 
Newtonian) limit. We will now derive a scalar-tensor 
representation for f(R) gravities that will allow us to 
treat the metric and Palatini formulations in a very sim- 
ilar manner and will simplify the computations of the 
post-Newtonian metric. The following action, leads to 
the same equations of motion as eq. © 

S = 2^ / <** V=5 if (A) + (R-A)B] + S m (2) 



where A and B are auxiliary fields. To fix our notation, 
we will denote R(g) the contraction R(g) = g^R^v with 
the Ricci tensor R^ v given in terms of the Levi-Civita 
connection of g^ v , and R(T) the contraction R(T) = 
g^Rfiv with R^y given in terms of a connection T in- 
dependent of g^y. Thus, the symbol R in eq.© must 
be seen as R = R(g) in the metric formalism, and as 
R = R(T) in the Palatini formalism. It is easy to see 
that the functional variation of eq. (J2J with respect to A 
gives B — df{A)/dA = 0. This equation can be used to 
algebraically solve for A — A(B). Inserting this solution 
back into eq.© and defining 



V(B)=AB-f(A) 



(3) 



then eq. J2J) in the metric formalism can be identified with 
the case oj — of the Brans-Dickc-likc theories 
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where <p = B, Analogously, one can find a scalar-tensor 
representation for the Palatini formulation. In this case, 
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the equations of motion for the connection lead to 



(5) 



where the tensor i M „ is defined as t^ v = Bg^. This 
solution for the connection allows us to write R(T) in 
terms of the metric and the field B as follows 

R{T) = R{g) + ^d x Bd x B-^UB (6) 

Inserting this solution in eq.@ and discarding a total 
divergence we get 
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which represents the case to = —3/2 of the theories de- 
fined in eq.|@}. It is worth noting that when V(B) is 
given, the inverse problem of finding the lagrangian f(R) 
is also possible. From eq.Q we see that 



dV{B) 
dB 
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Using this algebraic equation to solve for B — B(A), the 
lagrangian can be written, using again eq.©, as 



f(A) = AB — V(B) 



(9) 



This means that, with regard to the equations of motion, 
f(R) gravities in metric and Palatini formalisms are 
equivalent to ui = and lu = —3/2 Brans-Dickc-like 
theories respectively. 

The identification of f(R) gravities with particular 
cases of Brans-Dicke-like theories should, in principle, al- 
low us to express their post-Newtonian limit in a com- 
pact form dependent on the parameter lu. However, this 
is only partially true. A glance to the equation of motion 
for the scalar field defined in eq.Q 



[3 + 2w]D0 + 2V{(p) - (j>V'{<t>) = k2t 
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where T = g^T^, indicates that the field is a dynamical 
object for all u> but lu = —3/2. Thus, there exists a clear 
dynamical difference between a generic lu = constant 
theory and the case uj — —3/2, which corresponds to 
the Palatini form of f(R) gravities. As a scalar-tensor 
theory, the case uj = —3/2 seems to have been almost 
avoided in the literature. In the original Brans-Dicke 
theory, where the potential term was not present, this 
case was obviously pathological (see for a discussion 
of the limit uj — ► —3/2). Moreover, it was found that 
in order to get agreement between predictions and solar 
system experiments u> should be large and positive, 
and little attention was paid later, when non-trivial 
potentials were considered, to small or negative values 
of uj. on the other hand, for uj ^ —3/2 theories with 
V ^ it is well known that if the field has associated a 



large effective mass, the predictions of the theory may 
agree with solar system experiments irrespective of the 
value of w |||. However, such result assumes that the 
field is near an extremum of its potential. For uj = 
(f(R) in metric formalism) that condition, dV/d(f> = 0, 
cannot be imposed in general. This follows from the 
variation of eq.(J2J) with respect to B and eq.JHJ, which 
lead to dV/d(j) = R{g). Since the leading order of R(g) 
at a given time coincides with the background cosmic 
curvature Rq ^ 0, we cannot impose dV/d<p = for 
uj = theories. We are thus forced to compute the 
post-Newtonian limit of this case without making any a 
priori assumption or simplification about the behavior 
of the potential. In the case uj = —3/2 the relation 
between <\> and V(4>) (see ea. HlOJl ') is even stronger and 
also forces us not to make any assumption on V(<p). 

We shall now sketch the basic steps to compute 
the post-Newtonian metric of Brans-Dicke-like theories, 
which will be detailed elsewhere. We generalize the re- 
sults of the literature so as to include all the terms that 
are relevant for our discussion. For approximately static 
solutions, corresponding to masses such as the Sun or 
Earth, to lowest order, we can drop the terms involving 
time derivatives from the equations of motion. In a coor- 
dinate system in free fall with respect to the surrounding 
cosmological model, the metric can be expanded about 
its Minkowskian value as g^ v = rj^ u + h^ v . 
For uj ^ —3/2, the field is dynamical and can be ex- 
panded as 4> = 0o + up{t,x), where 4>q = 4>(to) is the 
asymptotic cosmic value, which is a slowly-varying func- 
tion of the cosmic time to, and ip(t, x) represents the local 
deviation from cf>Q. To lowest order (T ss — p), the metric 
satisfies the following equations 
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where the gauge condition h% — i/i^ k — duf^ /0o has 
been used. In eliminating the zeroth-order terms in the 
field equation for tp, corresponding to the cosmological 
solution for </>o, we obtain 
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where m 2 v is a slowly- varying function of the cosmic time 
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Solving eas. ((TT|l . (fT2"|) and (fH^ . the metric becomes 
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where Mq = J d 3 x'p(t, x') and we have defined the effec- 
tive Newton's constant and the PPN parameter 7 as 



G = 



F(r) 
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The function F(r) is given by F(r) 



(17) 
(18) 
when 



m 2 > and by F(r) 



cos |m v |r when < 0. 



Let us now obtain the metric for uj = —3/2 theories. In 
this case no boundary conditions are needed for the scalar 
field, since it satisfies an algebraic equation. Denoting 
by — 4>(T) the solution to ea. l|lU|) when u> = —3/2, 
we construct the quantity fi(T) = log[0/0o], where the 
subindex now denotes vacuum value, 0o = 4>{T = 0). 
Using the gauge condition h£ — k = 9/-Q(T), the 
equations for the metric can be written as 
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The solution to these equations are 
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where M Q = 4>o J d 3 x'p(t,x')/(j) and we have defined 

2 ' My\ 
M J 

Tl 77- ( 24 ) 



G 



87T0O 
Mq - My 



(23) 



Mq + My 

with My = K - 2 o /dV[y o /0 o - ^(0)M 

The term (Vb/0o)?" 2 appearing in eas. (|15|l and 116|) for 

w 7^ -3/2 and in and for w = -3/2 is related 
to the scalar energy density and acts in a manner simi- 
lar to the effects of a cosmological constant. Thus, any 
viable theory must give a negligible contribution of this 
type. Let us now focus on u> = theories. The oscillat- 
ing solutions F(r) — cos |m,Jr are always unphysical. If 
\m v \L <C 1 with L large compared to solar system scales 
(long-range interaction), we find 7 ss 1/2, which is ruled 
out by observations (jobs ~ 1, see 6]). If \m v \L > 1 
(short-range), Newton's constant strongly oscillates in 
space and the Newtonian limit is dramatically modified. 
Thus, only the damped case F(r) — e~\ m ^ T with 



m 2 L 2 > 1 



(25) 



is physically acceptable. This condition means that the 



scalar interaction range L 



must be shorter than 



any currently accessible experimental length L. 
In the case lo = —3/2, we can learn about the depen- 
dence of on T by studying the behavior of Mq , G and 
7. According to the definition of Mq and My (see below 
eas.l|22[l and Q24[0. it follows that a body with Newto- 
nian mass Mjv = / d 3 x'p(t, x') may yield different values 
of Mq, G and 7 depending on its internal structure and 
composition. In consequence, a given amount of Newto- 
nian mass could lead to gravitational fields of different 
strengths and dynamical properties. We thus demand a 
very weak dependence of on T so as to guarantee that 
Mq, G and 7 are almost constant. This is consistent with 
the requirement that the local term O(T) in eqs.JSTJ and 
1221 must be small compared to unity. Since the contri- 
bution of f2(T) to the acceleration of a body is given in 
terms of its gradient, we must demand that 



T{dcj>/dT) 



< 1 
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from T = up to nuclear densities, otherwise, a change 
in when going from its vacuum value 0o outside atoms 
to its value inside atoms would lead to observable effects 
in the motion of macroscopic test bodies placed in the 
gravitational field defined by eas. (|2"T]l and J32J. Since 
such effects have not been observed, it follows that over 
a wide range of densities w O + (d(j)/dT)\ T=0 T , with 
|0 o ~ 1 (90/9T)|t=oT| -C 1, must be a very good approxi- 
mation. The weak dependence on T expressed by ea- H26|) 
can be written using ea. 1)10(1 to evaluate (<90/<9T) as 
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We can thus interpret this equation in a manner analo- 
gous to eq. (|25|) , as the ratio of a length associated to the 



matter density, L 2 (p) — [n 2 p/(j)( 
ated to the scalar field, 2| 



, over a length associ- 

[(0y"-v")0/0o]- 1 . 



The constraint on u> = theories given in ea. l|25|l can 
be rewritten in terms of the lagrangian f(R) as follows 



f'(Ro) 
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where Rq represents the current cosmic scalar curvature. 
Note that since = /' > to have a well posed theory, 
it follows that /" must be small and positive in order to 
satisfy ea. l|28l) . In consequence, any theory with /" < 
leads to an ill defined post-Newtonian limit (m 2 < 0). 
This is the case, for instance, of the Carroll et al. model 
0- We shall now demand that the interaction range of 
the scalar field remains as short as today or decreases 
with time so as to avoid dramatic modifications of the 
gravitational dynamics in post-Newtonian systems with 
the cosmic expansion. This can be implemented imposing 
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as R — > 0, where P <§; £ 2 represents a bound to the cur- 
rent interaction range of the scalar field. Manipulating 
this inequality we obtain 



d\og[f(R)} 



< 



P 



dR ~ 1 + l 2 R 
which can be integrated twice to give 

l 2 R 2 



(30) 



f(R) <a + P[R 



where (3 > 0. Since /' and /" are positive, the lagrangian 
is also bounded from below, say, f(R) > a. According 
to the cosmological data, a = — 2A must be of order a 
cosmological constant 2A ~ 10~ 53 m 2 . Without loss of 
generality, setting (3 — 1 we find that in order to satisfy 
the current solar system constraints, the lagrangian in 
metric formalism must satisfy 



-2A < f(R) <R-2A 



l 2 R 2 



(32) 



which is clearly incompatible with nonlinear terms 
growing at low curvatures. 

Let us consider now the Palatini case, u> = —3/2. Writ- 
ten in terms of the lagrangian f(R), ea. (|27|l turns into 
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where /' = /'//o = 0/</V Since C 2 {p) ~ l/p takes its 
smallest value for ordinary matter at nuclear densities, it 
is reasonable to demand that 
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where I 2 represents a length scale much smaller than 
£ 2 (p) at those densities. Note that I 2 determines the 
scale over which the nonlinear corrections are relevant. If 
I 2 = 0, which implies f(R) = a + bR, then ea. l|33|l would 
be valid for all p. Furthermore, if the nonlinear correc- 
tions were important at very low cosmic densities, I would 
be of order the radius of the universe and the nonlinear 



terms would dominate at all scales, which would lead to 
unacceptable predictions, as we pointed out above with 
regard to M Q , G, 7 and fi(T). A good example of such 
pathological effects in the Palatini formalism was studied 
in [3 ■ Manipulating eq. l|3"4"|) we obtain for /" > 



f<a- 



l 2 R 2 R 
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(31) and for /" < 



pB 2 R 
— +2^ 



(l 2 Rf 



— \ og [l 2 R+^l + (PR) 2 ] 
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where the vacuum value f' has been set to unity and a 
can be identified with a = — 2A. We see that, to leading 
order in PR, the Palatini lagrangian is bounded by 

l2n2 PR 2 

R-2A--^<f(R)<R-2A+-^- (37) 

According to eqs.(j22l an d l|37|l . our conclusion is clear, 
laboratory and solar system experiments indicate that 
the gravity lagrangian is nearly linear in R, with the 
possible nonlinearities bounded by quadratic terms. In 
consequence, f(R) lagrangians with nonlinear terms that 
grow at low curvatures cannot represent a valid mecha- 
nism to justify the cosmic accelerated expansion rate. 
Such theories lead to a long-range scalar interaction in- 
compatible with the experimental tests. In the viable 
models the non-linearities represent a short-range scalar 
interaction, whose effect in the late-time cosmic dynam- 
ics reduces to that of a cosmological constant and, there- 
fore, do not substantially modify the description pro- 
vided by General Relativity. To conclude, we want to 
remark that eas. (|3^|l and (|33|l - l|S6^l tell us how the la- 
grangian f(R) must be near the origin, PR <C 1, not far 
from it, PR ^> 1, where the post-Newtonian constraints 
could not make sense. This fact constraints the possi- 
ble f{R) early-time inflationary models. In particular, 
f(R) = R + aR 2 in metric formalism seems compatible 
with CMBR observations [Toj . 
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